Towards an Understanding of Atmospheric Balance 


Ronald M. Errico 

Global Modeling and Assimilation Office 

GSFC NASA 

Goddard Earth Sciences Technology and Research Center 

Morgan State University 

(Formerly at) National Center for Atmospheric Research 



MORGAN STATE UNIVERSITY 



GESTAR 

v 


EXCELLENCE 


Why is the extra-tropical atmosphere approx, quasi-geostrophic 


The stability of quasi-geostrophic flow with respect to 

ageostrophic perturbations 


Derivation of (2-layer, f_plane) PE in terms of Normal Modes 
ErricoJAS 1981 
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Demonstartion of Equipartition 
Errico Tellus 1984 



Examination of Balance 
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NATURAL FREQUENCY |cr| (s' 1 ) 

Errico 1984, 1990 MWR: Errico et al. 1988 MWR 


Diabatic Balance ? 


The interplay of analysis and initialization 
Errico et al. MWR 1993 
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Gravitational modes considered as forced and damped harmonic oscillators 


Define 
time t , 
Then 


g(t) as the complex amplitude of a gravity- wave like mode at each 
and let B and G be the sets of Rossby- and gravity-wave like modes. 

^ = ~i\g + N(R) + N(R, G ) + N(G) + D(R, G) - vg 


Consider N(B) = F(t) as the dominant nonlinear term. Approximately then 


dg 

dt 


-i\g + F(t) - vg 


Consider F{t) = F(0) exp(— ipd). Then 
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Errico 1981 JAS, 1984 MWR, 1997 / Japan MS 


Harmonic Dial for External m=4 Mode, Period=3.7h 
Without NNMI With NNMI 




Errico 1997 J Japan Met Soc 


Behavior of gravitational modes in a climate model: 

Time series (harmonic dials) of complex mode amplitudes 


Errico MWR 1989 


16 days shown 
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Behavior of gravitational modes in a climate model 
Power spectra of complex mode amplitudes 
Errico MWR 1989 


m- 7 n= 4 i-1 WG p, = -3.1 
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Behavior of gravitational modes in a climate model 
Power spectra of convective heating 
Errico MWR 1989 
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Diabatic balance vs appropriate cutoff 
Errico and Rasch Tellus 1988 
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RELATIVE ACCUMULATED ERROR 


Higher-order Machenhauer schemes 
Errico MWR 1989 


ORDER 1 0-0 



ORDER 2 v- 0,0 
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Other Issues 


Vertical modes in discrete models 


10 level MAMS 72 level GEOS-5 

Modes 1, 2, 7 (H= 10,000, 2050, 13 m) Mode 29 (H=13m) 

(Notice only plotted up to o=0. 1 




23 zero crossings above for g<0. 1 


High amplitude modes in the upper atmosphere 


72 level GEOS-5 model with top at 0.01 hPa 


Vertical Modes: 0 to 3 



Global mean squared divergence tendency 
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Vertical Modes: 0 to 3 



Structures of 3 largest 
scale vertical normal 
modes 


GMAO-GSI 3DVAR 
72 level model 



Derivation of (2-layer) PE in terms of Normal Modes 
ErricoJAS 1981 
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FKt . 3. 1.ike Fig, l except for the fifth antisymmetric, zonal wave- 
number 4 P / - 5 WG mode. The progression i ti the dial is prcdonv 
inantly clockwise (west wand propagation). 
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Derivation of (2-layer) PE in terms of Normal Modes 
ErricoJAS 1981 
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Figure 4. The kinetic energy (K) T available potential energy (A), rotational -mode energy (R), gravitational -inode 
energy (G) + and total energy (E) contributed by vertical modes of indicated equivalent depths at t = 0_ The integers 
on the right-hand side indicate corresponds ng vertical-mode indices £ . 


Derivation of (2-layer) PE in terms of Normal Modes 
ErricoJAS 1981 




Figure 6. The (a) and (b) R and (c) and (d) its complement components of the (a) and (c) a' and (b) and <d) T' 
Helds on cr = 0.55 at t = 0 for SVl determined using the E norm applied to the dry form of the linearized model 
Contour intervals are (a) and (c) I ms' 1 , (b) I K, and Id) 0.5 K, with zero-contours omitted and negative values 

shown dashed. See text for further explanation. 


Derivation of (2-layer) PE in terms of Normal Modes 
ErricoJAS 1981 




Figure 8. The (a) and (c) t?' and (b) and! (d) T' fields on a = 0.55 at t = 24 h determined from the linearized 
evolutions begun from (a) and! (b) R-mode components of SV1 and (c) and (d) (heir complement of SVL Contour 
intervals are (a) 10 m s -1 , (b) 2 K f (c) 5 m s^ 1 . and (d) 1 K, with zero-contours omitted and negative values shown 

dashed. See text for further explanation. 



Partitioning of analysis error energy in terms of normal modes: 
(as inferred from an OSSE) 

Errico et al. Met Z 2007 
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Summary 


1 . Much can be learned from some old works 

2. The standard Normal Modes provide useful concepts and tools 

3. The standard Normal Modes have limitations 

a. the universality of vertical modes 

b. internal modes (when C approx = U) 

c. more realistic basic states (e.g. as for SVs) 

4. Is Initialization still an issue ? 

5. There is more to understand 

a. time scales of moist diabatic processes 

b. effects of top boundary conditions, non-hydrostatic behavior 

c. SY behavior 


